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Abstract 

This Letter reports evidence of triple gauge boson production pp —> W(/^v)j'y + X, which is 
accessible for the first time with the 8 TeV LHC data set. The fiducial cross section for this 
process is measured in a data sample corresponding to an integrated luminosity of 20.3 fb”L 
collected by the ATLAS detector in 2012. Events are selected using the W boson decay to 
ev or pv as well as requiring two isolated photons. The measured cross section is used to set 
limits on anomalous quartic gauge couplings in the high diphoton mass region. 


© 2015 CERN for the benefit of the ATLAS Collaboration. 

Reproduction of this article or parts of it is allowed as specihed in the CC-BY-3.0 license. 




In the Standard Model (SM), the self-couplings of the electroweak gauge bosons are specified by the non- 
Abelian S U(2) x U{1) structure of the electroweak sector. Since any deviation in the self-couplings from 
this expectation indicates the presence of new physics phenomena at unprobed energy scales, the meas¬ 
urement of the production of multiple electroweak gauge bosons represents an important test of the SM. 
This Letter presents a measurement of the triboson production cross section, discussed in Ref. [1], where 
the W boson decays into ev or fiv {W{€v)yy), and its sensitivity to anomalous quartic gauge couplings 
(AQGCs) VTVTyy. Such final sfafes mainly come from evenfs where fhe W boson is produced in fhe hard 
inferacfion befween fhe fwo parfons, and fhe phofons eifher originate from initial or final slate radiation 
processes, or from Iriple or quartic gauge vertices fogelher wilh fhe W boson. The inclusive and exclusive 
cross secfions are bolh measured. The inclusive case has no reslriclion on fhe Wyy recoil syslem, whereas 
fhe exclusive case includes a vefo on evenfs conlaining one or more jels. Limifs on AQGC parameters are 
sel in fhe exclusive phase space wifh a dipholon mass larger fhan 300 GeV. Tofal and differential cross 
secfions for fhe diboson produclion processes WW, WZ, ZZ, ITy, and Zy have been reported previously 
by fhe ATLAS [2-5], CMS [ 6 - 8 ], DO [9-11], and CDF [12-14] collaborations, including limifs on an¬ 
omalous Iriple gauge boson couplings. Limifs have been sel on AQGCs by ATLAS [15], CMS [16,17], 
fhe LEP experimenls [18-21], and DO [22]. 

ATLAS [23] is a mullipurpose defector composed of an inner fracking defeclor (ID) surrounded by a 
Ihin superconducling solenoid providing a 2 T axial magnetic field, elecfromagnefic (EM) and hadronic 
calorimelers, and a muon speclromeler (MS) immersed in fhe magnetic field produced by a syslem of 
superconducting loroids. Evenfs in Ibis analysis are selected wilh Iriggers requiring fhe presence of one 
muon wilh a Iransverse momenlum (px) of more fhan 18 GeV and fwo elecfromagnefic objecls wifh 
a Iransverse energy (Ex) of more fhan 10 GeV each, wilh an efficiency of aboul 80% [24], or Ihree 
Ex > 15 GeV elecfromagnefic objecls wilh an efficiency of more than 95% [25] [26]. After applying data 
quality requirements, the data set corresponds to a total integrated luminosity of 20.3 ± 0.6 fb~' [27]. 

The main backgrounds to the W{iv)yy process originate from processes with jets identified as phofons or 
leptons, referred lo as fakes hereafler. Dala-driven lechniques are used to eslimale fakes, whereas Monte 
Carlo (MC) simulafion is used fo eslimale background sources wilh prompl lepfons and photons and for 
the signal. The SHERPA 1.4.1 generator [28-31] is used to model the signal with up to three partons in 
the final slate. SHERPA was also used fo simulafe fhe Zy, Zyy, WZ, and lT(Tv)yy backgrounds. Eor fhe 
Zy background, fhe agreemenl befween dala and fhe MC prediction was assessed in Z-enriched conlrol 
regions. The tt, single top, and WW processes are modeled by MC@NEO 4.02 [32, 33], interfaced to 
HERWIG 6.520 [34] for parlon showering and fragmenlalion processes and fo JIMMY 4.30 [35] for 
underlying evenl simulafion. The POWHEG [36] generator is used lo simulate ZZ produclion, inlerfaced 
to PYTHIA 8.163 [37] for parlon showering and fragmenlalion. The CTIQ parlon disfribulion function 
(PDE) sel [38] is used for all SHERPA, MC@NLO, and POWHEG samples. The slandard ATEAS 
defeclor simulafion [39] based on GEANT4 [40] is used. If includes mulliple prolon-prolon inleraclions 
per bunch crossing (pile-up) as observed in dala. 

The W{iv)yy candidale evenfs conlain an isolafed lepfon and missing Iransverse momenlum (E™®®) from 
fhe undelecled neulrino of fhe leplonic W decay, and fwo isolafed phofons (including also phofons lhal 
have converted in electron-position pairs within the ID volume). Muon candidates are identified, wilhin 
pseudorapidify [26] \r]\ < 2.4, by associaling complefe tiacks or Irack segmenls in fhe MS wilh Iracks in 
fhe ID [41]. Election candidates are reconstiucfed wilhin \r]\ < 2.47 as elecfromagnefic clusters associ- 
aled to a Irack [42], whereas phofons are reconstiucfed as elecfromagnefic clusters wilh \r]\ < 2.37 [43]. 
The calorimeler Iransilion regions al 1.37 < I 77 I < 1.52 are excluded for elections and photons. Identi¬ 
fication criteria based on shower shapes in the EM calorimeter for photons, and additionally on tracking 
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Electron channel 

Muon channel 

Electron channel 

Muon channel 



> 0 

Ye, 

= 0 

Wj] + Wjj 

15.3 ± 4.8 (stat.) ± 5.3 (syst.) 

30.5 ± 7.7 (stat.) ± 6.8 (syst.) 

5.8 ± 2.1 (stat.) ± 2.0 (syst.) 

14.4 ± 4.9 (stat.) ± 4.9 (syst.) 

yy + jets 

1.5 ± 0.6 (stat.) ±1.0 (syst.) 

11.0 ± 4.0 (stat.) ± 4.9 (syst.) 

0.2 ± 0.2 (stat.) ± 0.2 (syst.) 

6.1 ± 3.5 (stat.) ± 3.1 (syst.) 

Zy 

11.2 ± 1.1 (stat.) 

3.9 ± 0.2 (stat.) 

2.4 ± 0.5 (stat.) 

2.8 ± 0.2 (stat.) 

Other backgrounds 

2.2 ± 0.6 (stat.) 

6.7 ± 2.0 (stat.) 

0.3 ±0.1 (stat.) 

1.1 ±0.3 (stat.) 

Total background 

30.2 ± 5.0 (stat.) ± 5.4 (syst.) 

52.1 ± 8.9 (stat.) ± 8.4 (syst.) 

8.7 ± 2.2 (stat.) ± 2.0 (syst.) 

24.4 ± 6.0 (stat.) ± 5.8 (syst.) 

Data 

47 

no 

15 

53 


Table 1; The background composition in each channel is shown for the inclusive (left) and exclusive (right) cases. 
The Wyi + VTjj and yy + jets backgrounds are estimated using data-driven techniques, whereas the others are 
extracted from MC simulation. The number of candidate events in data passing the full selection is also shown. 


information for electrons, referred to as “tight” in Refs. [42,44], are used. The uses the energy 
deposits in the calorimeters within \t]\ <4.9 and the muons identified in the MS, as described in Ref. [45]. 
Reconstructed muons, electrons, and photons are required to have > 20 GeV and to be isolated. 
Photons are considered isolated if the sum of calorimeter transverse energy deposits in a cone of size 
AR = 0.4 around the candidate is smaller than 4 GeV. The isolation is corrected for photon energy leak¬ 
age. The muon isolation is based on the sum of the transverse momenta of ID tracks in a cone of size 
AR = 0.2 which must be below 0.15 x pj. For electrons, the calorimeter transverse energy deposits and 
the sum of the transverse momenta of tracks in a cone of size A/? = 0.2 must be below 0.2 x and 
0.15 X p^, respectively. The lepton must also be compatible with originating from the primary vertex of 
the interaction, which is taken to be the vertex with the largest of associated tracks. F"™®® is required 
to exceed 25 GeV. The transverse mass of the W boson [46] is required to be greater than 40 GeV. The 
two photons must be outside of their mutual isolation cones by requiring AR{y, y) > 0.4. To suppress the 
contribution from final-state radiation, the lepton and photons are required to have AR(^, y) > 0.7. Events 
containing a second reconstructed lepton are rejected to reduce background from Drell-Yan events. In the 
electron channel, additional requirements are used to suppress events in which one electron is misidenti- 
fied as a photon (mainly originated from the Zy process): the transverse momentum of the eyy system is 
required to be greater than 30 GeV, and the invariant mass of the electron and the leading, subleading or 
both photons is required to be outside a 13, 8 or 15 GeV wide window around the Z boson mass, respect¬ 
ively. Exclusive events are defined wifh a vefo on addifional jefs compared fo fhe inclusive selecfion. Jefs 
are reconsfrucfed from clusfered energy deposifs in fhe calorimeter using fhe anfi-^^ algorifhm [47] wifh 
radius parameter R - 0.4 and are required fo have px > 30 GeV and |? 7 | < 4.4. Jefs af A/? < 0.3 from fhe 
selected lepfon and photons are rejecfed. In order fo reduce pile-up effecfs, for jefs wifh pj < 50 GeV 
and \r]\ < 2.4, more fhan 50% of fhe summed scalar px of fracks wifhin AR = 0.4 of fhe jef axis musf be 
from fracks associafed fo fhe primary verfex. 

Table 1 shows fhe expecfed background as well as fhe observafion. The background expecfafion alone is 
nol sufficienf fo describe fhe dafa indicafing fhe presence of signal evenfs. The fake-photon background 
from VFyj -i- VFjj is esfimafed by performing a fwo-dimensional femplafe fil fo fhe isolation energy disfri- 
bufions of fhe leading and subleading photons, as described in Ref. [48]. Three background femplafes are 
obfained from dafa by reversing some of fhe phofon identification requiremenfs based on shower shape; 
fhe signal femplafes are faken from MC simulafion. Confribufions from evenfs where a jef satisfies fhe 
elecfron idenfificafion criteria, or fhe muon originafes from heavy-flavor decays, i.e. from yy -i- jefs pro¬ 
cesses, are esfimafed by using a fwo-dimensional sideband mefhod consfrucfed from fhe lepton isolation 
and variables, as described in Ref. [5]. The distribufion of fhe diphofon invarianf mass in fhe fwo 
channels is shown in Eig. 1. Alternative mefhods have been used fo cross-check fhe esfimafe of fhe back- 


3 













• Data 
I I Wyy 

I I Wyj + Wjj 

I I yy + jets 

Other backgrounds 

= 8 TeV, 20.3 fb ' 



electron channel (N > 0) 

' jets ' 


400 500 

[GeV] 



Figure 1; Diphoton invariant mass distribution in the electron (left) and muon (right) channels. The expected signal 
based on the SHERPA prediction is shown. The hashed areas show the total systematic and statistical uncertainty 
on the background estimate. 


grounds coming from fakes, all of them provide consistent results. In the estimation of the fake-photon 
background, systematic uncertainties arise from the limited number of events in the control regions, the 
functional form used to describe the background isolation energy distribution, the definition of the control 
region, the modeling of the signal in the MC samples and the corresponding statistical uncertainty. In the 
estimate of the fake-lepton background, systematic uncertainties related to the control region definitions 
and the residual correlation of the discriminating variables are considered. 

The fiducial cross sections are obtained from a maximum-likelihood fit, similarly to Ref. [5] , for the 
electron channel, the muon channel, and the combination of the two assuming lepton universality to de¬ 
termine the 'W{lv)yy cross section for a single lepton flavor. They are measured in a phase space, defined 
in Table 2, close to that of the experimentally selected region. Here pj is the transverse momentum of 
the neutrino and is the fractional energy carried by the closest particle-level jet in a cone of A/? = 0.4 
around each photon direction. 


Definition of the fiducial region 
> 20 GeV, p\ > 25 GeV, \rie\ < 2.5 
> 40 GeV 

> 20 GeV, \if\ < 2.37, iso. fraction < 0.5 
AR(Ay) > 0.7, AR(y,y) > 0.4, AR(7’/y,jet) > 0.3 

Exclusive: no anti-A:, jets with > 30 GeV, < 4.4 


Table 2: Definition of the fiducial region for which the cross section is evaluated. 


The efficiency of the signal selection and the small acceptance correction due to the extrapolation over 
the calorimeter transition region and to Ipl = 2.5 for the leptons are taken into account in the procedure. 
The acceptance correction factors are 0.83 and 0.90 in the electron and muon channel, respectively. The 
combined efficiency and acceptance correction amounts to (19.6 ±0.5)% and (40.4 ±0.7)% in the electron 
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[fb] 

^MCrM 

Inclusive (Ajet > 0) 


pvyy 

evyy 

ivyy 

7.1 ];j '2 (stat.) ±1.5 (syst.) ±0.2 (lumi.) 
4.3 g (stat.) g (syst.) ±0.2 (lumi.) 

6.1 (stat.) ±1.2 (syst.) ±0.2 (lumi.) 

2.90 ±0.16 

Exclusive (Ajet = 0) 


pvyy 

evyy 

ivyy 

3.5 ± 0.9 (stat.) ];} Q (syst.) ±0.1 (lumi.) 

1.9 +} | (stat.) i (syst.) ±0.1 (lumi.) 

2.9 (stat.) +1'° (syst.) ±0.1 (lumi.) 

1.88 ±0.20 


Table 3; Measurement of the pp ivyy + X inclusive and exclusive fiducial cross sections. 


and muon channels in the inclusive case, and to (15.1 ±0.7)% and (39.7 ± 1.0)% in the exclusive case. The 
given uncertainties are statistical only. Corrections are applied to account for small differences between 
data and MC simulation in lepton, photon, and jet efficiencies, momentum scale and resolution, additional 
pp interactions, and beam-spot position. 

Systematic uncertainties on the cross section are accounted for by introducing nuisance parameters in 
the likelihood which modify the signal and background expected yields. Correlations between systematic 
uncertainties in the two channels are accounted for in the combined fit. When combining the two channels, 
the dominant systematic uncertainties in the inclusive and exclusive cross-section measurements are 14% 
and 23% from the data-driven background estimates, 5% to 7% from the jet energy scale, and 3% from the 
luminosity. Other systematic uncertainties considered stem from the electromagnetic and muonic energy 
scale and resolution, the object reconstruction, the pileup description, and the trigger efficiency. These 
are found to have a minor impact, below 3%. Theoretical uncertainties on the signal modeling, affecting 
only the acceptance extrapolation, are negligible. The measured cross sections are shown in Table 3. The 
significance after combining the two channels is larger than 3cr in the inclusive case. The measurements 
in the electron and muon channels are compatible within Icr. 

The SM prediction for the W(£v)yy cross section is calculated with the parton-level MC program MCFM [49] 
at next-to-leading order (NLO). The calculations are performed using the MCFM default electroweak 
parameters [50] and the CTIQ PDF set. The renormalization and factorization scales are set to the in¬ 
variant mass of the ivyy system. The fragmentation of quarks and gluons to photons is included using 
the fragmentation function GdRG_L0 [51]. The kinematic requirements at parton level match the fiducial 
acceptance of Table 2. 

In addition to the inclusive prediction, an exclusive cross section is obtained by vetoing events with 
an additional jet emission. To account for the difference between jets defined at parton and particle 
levels, a correction factor of about 0.87 in the exclusive case is computed and applied to the prediction as 
documented in Ref. [5]. Uncertainties on the two predictions include the effect of varying independently 
the renormalization and factorization scales by factors of 0.5 and 2.0, evaluating the CTIQ PDF error sets 
scaled to the 68% confidence level (CL), the uncertainties on quark or gluon fragmentation to a photon, 
and the parton to particle correction factors. The predictions for W(iv)yy production are compared to 
the measured cross sections in Table 3. The measured cross section is higher by 1.9cr in the inclusive 
case, while better agreement is seen in the exclusive case, similar to the measurement of Wy and Zy 
in Ref. [5]. In the case of Zy and Wy, higher order corrections were calculated to be smaller for the 
exclusive compared to the inclusive case [52]. As the process Wyy has similar properties, the exclusive 
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measurement is expeeted to be in better agreement with the theoretical prediction than the inclusive one. 
Therefore, in the following, the exclusive measurement will be used for the AQGC limits setting, as done 
in Ref. [5]. 

The AQGCs are introduced as dimension-8 operators following the formalism defined in the Appendix 
of Ref. [53]. While many operators give rise to anomalous couplings of the form WVTyy, this study is 
restricted to /to/A^, /m 2 /A.^> and /ms/A"^, where A represents the scale at which new physics appears, 
and / the coupling of the respective operator. The Wyy final state is expected to be particularly sensitive 
to the TO operator, whereas the other two operators can be related to the parameters of the dimension- 
6 operators used at LEP [18-21] and by CMS [16] via the transformations described in Ref. [54]. To 
preserve unitarity up to high energy scales, a form factor is introduced which depends on the energy, the 
form factor scale Aff and an exponent n, following the formalism described in Refs. [55,56]. The scale 
Aff is independent of the new physics scale A [57]. The largest form factor scale ensuring unitarity for 
this process at = 8 TeV, calculated using the VBFNLO generator [58-61], is given by n = 2 and 
Aff = 600 GeV for and Aff = 500 GeV for and /ms/A"^. 

Deviations from the SM prediction for the AQGC parameters, which are predicted to be zero, lead to an 
excess of events with high diphoton invariant mass. The phase space to study AQGCs was optimized 
using the expected significance calculated on simulated events. The optimal phase space was found to 
be the exclusive selection with the additional requirement of tUyy > 300 GeV. The SM backgrounds 
in this region are determined from a fit to the observed ruyy distribution. The expected SM background 
is 0.01 ± 0.03 (stat.) ± 0.20 (syst.) ( 0.02 ± 0.05 (stat.) ± 0.46 (syst.) ) events in the electron (muon) 
channel, where uncertainties include systematic effects due to the extrapolation procedure, depending on 
the modelling of the spectrum and on the initial background estimate used in the extrapolation. No events 
are observed in the high-mass region. 




Observed [ TeV ^] 

Expected [ TeV ^] 

n = 0 

/to/A"^ 

[-0.9, 0.9] X 10^ 
[-0.8, 0.8] X 10^ 
[-1.5, 1.4] X 10"^ 

[-1.2, 1.2] X 10^ 
[-1.1, 1.1] X 10^ 
[-1.9, 1.8] X 10"^ 

n = 1 

/to/A'* 
/m2/A^ 
/m3/A^ 

[-7.6, 7.3] X 10^ 
[-4.4, 4.6] X 10^ 
[-8.9, 8.0] X 10^ 

[-9.6, 9.5] X 10^ 
[-5.7, 5.9] X 10^ 
[-11.0, 10.0] X 10^ 

n - 2 

/to/A^ 

/m2/A^ 

/m3/A^ 

[-2.7, 2.6] X 10^ 
[-1.3, 1.3] X 10^ 
[-2.9, 2.5] X 10^ 

[-3.5, 3.4] X 10^ 
[-1.6, 1.7] X 10^ 
[-3.7, 3.3] X 10^ 


Table 4: Observed and expected 95% CL limits obtained for the /to/A'*, /m 2 /A'* and /mb/A"^ AQGC parameters for 
the combination of the two channels. The values of n = 0, 1, 2 are the exponential choices of the form factor, App 
is fixed to 600 GeV for /to/V^ and to 500 GeV for the other parameters. The n - 0 choice produces the limits 
without the form factor applied. 

The cross-section prediction as a quadratic function of the AQGC parameters is obtained by using VB¬ 
FNLO [58-61]. For SM couplings VBFNLO agrees with MCFM. The limits on the AQGC parameters 
are extracted with a frequentist profile likelihood test [62], using the methodology of Ref. [5]. The ex¬ 
pected and observed limits at 95% CL on the AQGC parameters are shown in Table 4 for different values 
of n. The limits on /m 2 /A^ and /ms/A^ improve on the previous results from LEP [18-21] and DO [22], 
but are less stringent than those from CMS [16, 17]. The limit on /to /A"* is tighter than the previous 
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limit published by CMS [17,63]. This can be explained by the fact that /to/A^ is especially sensitive to 
transversely polarized W bosons, which are favored in the present study [53]. 

In summary, evidence for the W{iv)yy process is reported for the first time. The significance of fhe 
inclusive production cross section is larger fhan 3 cr. The measured cross sections are in agreemenf wifhin 
uncerfainfies wifh NLO SM predicfions calculafed wifh MCFM. Limifs are sef af 95% CL on fhe AQGC 
parameters, in particular improving fhe limif on /to/A^. 
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